Resonant Cooper Pair Tunneling: Quantum Noise and Measurement Characteristics 
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We study the quantum charge noise and measurement properties of the double Cooper pair reso- 
nance point in a superconducting single-electron transistor (SSET) coupled to a Josephson charge 
qubit. Using a density matrix approach for the coupled system, we obtain a full description of the 
measurement back-action; for weak coupling, this is used to extract the quantum charge noise. Un- 
like the case of a non-superconducting SET, the back-action here can induce population inversion 
in the qubit. We find that the Cooper pair resonance process allows for a much better measurement 
than a similar non-superconducting SET, and can approach the quantum limit of efficiency. 

PACS numbers: 73.20.Dx., 73.23.Hk, 73.40. Gk 
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Among the many open issues related to solid state 
quantum computation, the question of how best to mea- 
sure a solid state qubit remains a particularly interesting 
one. In the case where the qubit is a Cooper pair box (i.e. 
a Josephson-junction single charge box), the standard 
choice for a read-out device is the single-electron tran- 
sistor (SET) An alternate and potentially more 
powerful approach is to use a superconducting single elec- 
tron transistor (SSET) biased at a point where the cyclic 
resonant tunneling of Cooper pairs dominates transport 
[@-|llJ Such processes, known as Josephson quasi- 
particle (JQP) resonances, would appear to be an attrac- 
tive choice for use in a measurement as their resonance 
structure implies an extremely high sensitivity. However, 
precisely because of their large gain, these processes may 
be expected to strongly alter the state of the qubit in a 
measurement. To assess the balance between these two 
opposing tendencies, a close examination of the physics 
of JQP tunneling is required. Note that the SSET-qubit 
system is more relevant to experiment than a setup hav- 
ing a SET, as fabrication typically results in the entire 
system being made from a single metal. 

In this paper, we focus on a double JQP process 
(DJQP) (see Fig. g), which occurs at a lower SSET 
source-drain voltage than single JQP processes, and 
which has been used in a recent experiment |i~3| ]. We 
assess the potential of DJQP to act as a one-shot mea- 
surement of the state of a Cooper pair box qubit. This 
involves characterizing both r moas , the time needed to 
discriminate the two qubit states in the measurement, 
and the back-action of the measurement on the qubit, 
which is described by a mixing rate r m i X and a dephas- 
ing rate t v . These quantities are intimately related to 
the noise properties of the SSET, which are of funda- 
mental interest in themselves, given the novel nature of 
the DJQP process. r mcas is determined largely by the 
shot noise of the process, while r m i x and t v are related 
to the associated charge noise on the SSET island. While 
the shot noise of a single JQP process has been analyzed 
recently |lj] , the quantum charge noise has not been ad- 
dressed. It is of particular interest, as the experiment of 
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FIG. 1. Schematic showing the four steps of the double 
Josephson quasiparticle process which can occur in a super- 
conducting single-electron transistor. Circles represent the 
central island of the SSET, while the rectangles are the elec- 
trodes. Numbers indicate the charge of the SSET island. 

Rcf. jl3| uses the ability of a qubit to act as a spectrum 
analyzer of this quantum charge noise. 

To describe the measurement process in our system, 
we employ a density matrix description of the fully cou- 
pled SSET plus qubit system; this is similar to the ap- 
proach taken by Makhlin et al. J| for a SET, but ex- 
tended to deal with Josephson tunneling. This approach 
is not limited by a requirement of weak-coupling, as are 
standard approaches which perturbatively link r m ix to 
the transistor charge noise nonetheless, in the limit 
of weak-coupling the present method can be used to cal- 
culate the quantum charge noise of the SSET. We find 
that the quantum (i.e. asymmetric in frequency) nature 
of the noise is particularly pronounced for the DJQP fea- 
ture, leading to regimes where the SSET can strongly 
relax the qubit. Moreover, due to the resonant nature 
of Cooper pair tunneling, there exist regimes where the 
SSET can cause a pronounced population inversion in 
the Cooper pair box, something that is impossible us- 
ing a SET. For typical device parameters, we find that a 
far better single-shot measurement is possible using the 
DJQP process than with a comparable SET. 

Model- The Hamiltonian of the coupled qubit plus 
SSET system is is written as H — Hs + + TLmt- The 
qubit itself (or "box"), described by TLq, consists of a 
superconducting metal island in the Coulomb blockade 
regime where only two charge states are relevant. These 
can be regarded as the o~ z eigenstates of a fictitious spin 
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1/2. The island is attached via a tunnel junction to a 
bulk superconducting electrode, leading to the form 

n Q = ~ [{AE C q{1 - Nq)) a z + E JQ o x ] . (1) 

where Ecq is the charging energy of the box, Ejq is 
the Josephson coupling energy of the box, and Aq is 
the dimensionless gate voltage applied to the box. The 
SSET consists of a superconducting, Coulomb-blockaded 
island which is attached via tunnel junctions to two su- 
perconducting electrodes (Fig. |l|). The SSET Hamilto- 
nian Hs = Hk+Hc+Hv+Ht has a term Hk describing 
the kinetic energy of source, drain and central island elec- 
trons, a term Hy which describes the work done by the 
voltage sources, and a tunneling term Ht- The charging 
term is H c = Ecs(ns ~ As) 2 , where E C s is the SSET 
charging energy, n$ is the number of electrons on the 
central island, and As is the dimensionless gate voltage 
applied to the island. Finally, the qubit is capacitively 
coupled to the SSET: Hint = 2EcQ%£(?zns = E int a z ns- 
Here Cc is the cross-capacitance between the box and the 
central island of the SSET, and Cs is the total capaci- 
tance of the SSET island. Note that we neglect the cou- 
pling of the qubit to its environment, as we are interested 
here in the intrinsic effect of the SSET on the qubit p5| . 
We also assume a SSET with identical tunnel junctions, 
whose dimensionless conductance g satisfies g/(2ir) -C 1. 

The DJQP process occurs when the SSET gate volt- 
age As and drain-source voltage 2Vds are such that two 
Cooper-pair tunneling transitions (one in each junction) 
are resonant. We label these transitions as n s — — > 2 
(left junction) and n s = 1 — ► — 1 (right junction) (see 
Fig. 0). Resonance thus requires eVjjs — Ecs and 
As = 1/2. In addition, Ecs /As (where As is the super- 
conducting gap of the SSET) must be chosen so that the 
quasiparticle transitions linking the two Cooper pair res- 
onances are energetically allowed (i.e. ns — 2 — > 1 and 
ns = —1 — ► 0), whereas transitions which end the cycle 
(i.e. ns = — > 1) are not. We take Ecs — As to sat- 
isfy these conditions; this corresponds to the experiment 
of Ref. p]J . The two quasiparticle transitions which oc- 
cur in the DJQP are characterized by a rate T, which is 
given by the usual expression for quasiparticle tunneling 
between two superconductors ITS). The effective Cooper 
pair tunneling rate jj emerging from our description (i.e. 
Eq. (|J) below) is given by B: 



7J(5) 



4(<52 + (r/2)2) 



(2) 



Here, 6 is the energy difference between the two charge 
states involved in tunneling, Ejs is the Josephson energy 
of the SSET, and we set % = 1. 

Calculation Approach- We consider the reduced den- 
sity matrix p of the qubit plus SSET system obtained by 
tracing out the SSET fermionic degrees of freedom. The 



evolution of p is calculated perturbatively in the tunnel- 
ing Hamiltonian Ht, keeping only the lowest order terms; 
this corresponds to the neglect of co-tunneling processes, 
which is valid for small g and near the DJQP resonance. 
Using an interaction representation where TLt is viewed 
as a perturbation, the equation of motion of p takes the 
standard form: 



dt 



p(t) 



dt'{[H T (t),[HTtf),p(t')® PF ]]) (3) 



The angular brackets denote the trace over SSET fermion 
degrees of freedom; as we work at zero temperature, pf 
is the density matrix corresponding to ground state of 
these degrees of freedom in the absence of tunneling. In 
the diagrammatic language of Ref. , Eq. (||) is equiv- 
alent to keeping all Tt T terms in the self-energy of the 
Keldysh propagator governing the evolution of p. Note 
that the correlators in Eq. (||) describe both quasiparticle 
tunneling and Josephson tunneling in the SSET. 

To make further progress, we treat the Josephson cou- 
pling emerging from Eq. (^|) as energy-independent and 
given by the Ambegaokar-Baratoff value Ejs = gAs/8. 
We also use the smallness of g to neglect logarithmic 
renormalization terms, as was done in Ref. [0. One can 
then immediately solve for the time-independent solution 
of Eq. (||), which describes the quasi-equilibrium state 
achieved by the system after all mixing and dephasing 
of the qubit by the SSET has occurred. To describe the 
dynamics of mixing (i.e. the relaxation of the qubit state 
populations to their stationary value), we also calculate 
the corresponding eigenmode of Eq. (||) . A Markov ap- 
proximation is made which involves replacing p(t') by 
p{t) on the RHS of Eq. (|); for the mixing mode, this 
should be done in the Schrodinger picture |Q . This ap- 
proximation is justified as long as the time dependence 
of p in the mixing mode is weak compared to typical fre- 
quencies appearing in the correlators of Eq. (p[), requiring 
here that T mix < T, E C s and E JS < E C s @- 

Back-Action- We focus here primarily on the mixing 
effect of the measurement back-action; dephasing will be 
discussed more extensively in Ref. | fL8[ . The mixing rate 
Tmix = r rc i + r cxc is set by the rates at which the mea- 
surement relaxes and excites the qubit. Let denote the 
A/Q-dependent energy difference between the two qubit 
states. For weak coupling (Ei nt <C fi), Fermi's Golden 
rule relates r re i and r eX c to the quantum charge noise of 
the SSET island Sq{lj) = J dt e~ iujt (n s (t)n s (0)): 



r 



rcl/cxc 



= El 



5 Q (±n). 



(4) 



In our approach, these rates may be directly obtained 
by using the stationary solution (which gives the post- 
mixing occupancies of the box eigenstates) and the mix- 
ing eigenvalue of Eq. (3). In the limit of weak coupling, 
one can then use Eq. (4) to extract Sq(Q). Our method 
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FIG. 2. Quantum charge noise associated with the DJQP 
process. The solid curve corresponds to Ms, Vds tuned to the 
center of the DJQP resonance; the dashed curve corresponds 
to moving eVos away from resonance by +F/4. We take 
g — 0.5 and As ^ Ecs — 0.25meV in the SSET, correspond- 
ing to the device of Ref. 13; this gives Ejs / {hT) ~ 0.04. Inset: 
average qubit charge after mixing has occurred for weak cou- 
pling (Eint/EjQ = 0.01), as a function of qubit gate voltage 
A/q; see text for details. As in Ref. 13, we have Ecq — 77{ieV 
and Ejq ~ 27fieV. The frequency range probed by tuning 
Mq matches the range of the main plot; the sharp steps in 
the average charge occur at Q(Mq) — Ecs- 

for calculating the quantum noise, which uses the qubit as 
a spectrum analyzer, is physically intuitive and no more 
difficult to implement than standard approaches ||; in 
addition, we are able to calculate T re i and r exc when the 
coupling is not weak, and Eq. (||) fails. 

Fig. ^| displays the quantum charge noise obtained at 
zero temperature, using SSET parameters which corre- 
spond to Ref. [13]. The solid curve in Fig. ||| is for the 
center of the DJQP resonance- Ms = 1/2, eVbs = Ecs- 
Note the sudden asymmetry that develops between ab- 
sorption (i.e. Sq(-|-|w|)) and emission (i.e Sq(— \uj\)) 
when |cj| increases beyond Ecs- These sudden jumps 
correspond to the opening and closing of transport chan- 
nels in the SSET, and their sharpness (which is absent 
for similar processes in a SET Q) is a direct consequence 
of the singularity in the quasiparticle density of states. 
For example, as uj rises past EcSi quasiparticle transi- 
tions which are normally forbidden in the DJQP cycle 
(i.e. ns — — » 1) suddenly become energetically allowed 
if they absorb energy from the qubit, causing a sudden 
increase in Sq(uj). 

The effect of the SSET quantum charge noise on the 
qubit is shown in the inset of Fig. |^, where the average 
qubit charge (Ng) = 1 + (cr 2 ) for t 3> T mix is shown 
as a function of A/q. Changing A/q tunes the qubit 
splitting frequency tt, allowing one to effectively probe 
the frequency dependence of the noise. The solid black 
curve corresponds to being at the center of the DJQP 
feature, while the grey curve corresponds to the unper- 
turbed qubit ground state. The sharp features in the 
quantum noise manifest themselves in (Ng), a quantity 
which is directly accessible in experiment. 

Even more interesting are features emerging in the low 
frequency quantum noise (|tt>| <C Ecs) when one tunes 



Ms or Vds slightly off the DJQP resonance center. Un- 
like the case of a SET, where asymmetries in the noise 
are weak for these frequencies, there are strong features 
here that result from the resonant nature of Cooper pair 
tunneling. By treating the mixing terms in Eq. (^) per- 
turbatively, simple analytic expressions can be obtained 
for the quantum noise in this regime when Ejs < T (in 
Ref. H|, Ej S /(hT) ~ 0.04). If one moves away from 
the DJQP center by tuning only Vds (i-e. Ms = 1/2, 
eV DS = E C s + <V/2), we find (|w| < E C s)- 



[u(o~v + u)hj(6y - ui)} 
[4r/j(6 v + Lu)-f,j(5 v - uj)] + uj 2 



(5) 



In the limit where uj is much smaller than the width T/2 
of the Cooper pair resonance, Eq. (||) simply corresponds 
to classical telegraph noise (the SSET only spends appre- 
ciable time in the states ns = and ns = 1). However, 
for finite Sy and uj, Eq. (||) indicates that the noise de- 
velops a pronounced asymmetry, even though \uj\ <§C Ec- 
In particular, if Sy > 0, one has Sq(—\uj\) > Sq(+\uj\), 
implying that emission by the SSET exceeds absorption. 
This behavior is shown by the dashed curves in Fig. |^, 
which correspond to Ms = 1/2, S v = +T/4. This ef- 
fect is a direct consequence of the resonant nature of 
Cooper pair tunneling- by emitting energy, both Cooper 
pair tunneling processes in the DJQP cycle become more 
resonant, while absorbing energy pushes them even fur- 
ther from resonance. The net result is a population in- 
version in the qubit at zero temperature, which in turn 
leads to a striking, non-monotonic dependence of qubit 
charge on Mn (this is shown by the dashed curve in the 
inset of Fig. g) jTJ) . Note that if one moves away from 
the center of the DJQP resonance by only changing the 
gate voltage Ms, no asymmetry in the noise results, as 
now emission (or absorption) moves one of the Cooper- 
pair transitions in the DJQP process further towards res- 
onance, while it moves the other transition further away 
from resonance. Nonetheless, the noise in this case still 
has a non-monotonic dependence on frequency. Letting 
S v = and 5jy = AE C s{Ms - 1/2), we have for E JS < T: 

l+(Mpjt 7J ( S/V -uj) 7 j(5M + uj) 

S q (uj)= 7 j(S /V ) jf- ^— — (6) 

[47j((W + uJ)^j(dM - uj)\ + uj 

Measurement Rate- To determine the measurement 
time Tmeas, we extend our density matrix description to 
also include m, the number of electrons that have tun- 
nelled through the left SSET junction @,@ ||. We are 
thus able to calculate the distribution of tunnelled elec- 
trons P(m,t\i), where i =T, J, denotes the initial state of 
the qubit. r mcas is defined as the minimum time needed 
before the two distributions P(m, t\ |) and P(m, t\ 1) are 
statistically distinguishable [I): 



72/Vt + VWT 



(J) 
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Here, 7f and I± are the average SSET currents associ- 
ated with the two qubit states, and /j and f± are the 
associated Fano factors which govern the zero-frequency 
shot noise in the current (i.e. the distribution P(m,t\i) is 
roughly Gaussian with a mean (Jije)t and standard de- 
viation J fi(Ii/e)t). In the absence of the qubit, analysis 
of the density matrix equations for the SSET reveals: 



f(s) = - 2 



E JS 



1 

1_ 2 ([ r/2p 



(3(r/2) 2 - s 2 



s 2 



(8) 



where we take eVos = E C s, 5 — <W = AE C s{Ns - 1/2). 
Eq. (g) indicates that the effective charge of the carriers 
in the DJQP process is 3e/2 in the limit where Y ^> Ejs- 
In this limit, Cooper-pair tunneling is the rate-limiting 
step in the cycle; electrons effectively tunnel in clumps of 
e or 2e, leading to an average charge of 3e/2. A similar 
argument holds in the opposite regime Y <C Ejg- 

We consider r moas in the limit of weak coupling (E int <C 
Yl) and weak mixing (Ejq <C Yl). In this limit, the two 
qubit states simply provide an effective shift in the SSET 
gate voltage. Taking 8y — and — Y/2 for near opti- 
mal gain, and using Eqs. (f|Mh we find that the intrinsic 



signal-to-noise ratio (r meas r D 



-1/2 



of the measurement, 



in the relevant regime Ejs < Y, is given by: 
(S/N) D j Q p = ^\cot9\^. 



(9) 



Here, cot 6 = <±E C q(1 -Nq)/E jq , and we take jj(0) < 
Yl < Ecs- If a SET in the sequential-tunneling regime 
is used for the qubit measurement, it was found in Rcfs. 
that the optimal S/N is given by (Yl < Ecs)'- 



(S/N)set — A I cot t? I 



eVos 
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ZS, (10) 



where A is of order unity. As the quasiparticle transition 
rate Y ~ j-eVoS; we see that the S/N achieved using 
DJQP is paramctrically larger (in 2ir/g ^> 1) than that 
obtained for the SET. This enhancement results largely 
from the narrow width of the DJQP feature- the energy 
scale over which the current changes (and thus the gain) 
is set by Y rather than Vds- The gain and S/N ra- 
tio of the SET could be improved by working in the co- 
tunneling regime; however, this would result in a much 
larger T rneas (r rneas oc g~ 2 ), making one more suscepti- 
ble to unwanted environmental effects. In contrast, the 
DJQP feature has both a large gain and a short T meas 
(i.e. T meas oc 1/g). Shown in Fig. | as a function of 
Nq are T mea s,rrci and r oxc for a strongly coupled device 
(Cc/Cs = 0.05), with all other parameters correspond- 
ing to Ref. E|. We have taken Sy = and 6tf = Y/2 for 
optimal gain. Fig. ^ confirms that an excellent measure- 
ment is indeed possible, with (S/N) 2 > 100. Note that 
the strong coupling splits the sharp features occurring in 
the noise (see Fig [| solid curve). 




Qubit Gate Voltage 2 

FIG. 3. l/r mcas , r rc i, and r cxc vs. qubit gate voltage 
A/q for a strongly coupled system, where Ei nt /EjQ ~ 0.3 
(i.e. Cc/Ce = 0.05). A good measurement is possible for 
a wide range of gate voltages. Inset: Heisenberg efficiency 
X = 7^/rmeas at weak coupling, as a function of Ejs/Y. 

We have also studied the "Heisenberg efficiency" \ — 
T v /r roea ,s of measurement using DJQP for a weak cou- 
pling (E'int < E JS ,Y) and Yl < E C s, where t v is the 
measurement-induced dephasing time [ jl8| . Unlike an 
SET in the sequential tunneling regime, where x x 9 2 
is always much less than the quantum limit x — 1 
here x 1S controlled by the ratio Ejs/Y. As shown in the 
inset of Fig. |^, by tuning this ratio, x can be made to ap- 
proach the quantum limit. Here, for each value of Ejs /Y 
we have set Vds an d N"s to optimize the gain. Measure- 
ment using DJQP is able to reach a high efficiency when 
Ejs — Y both because of the symmetry of the process, 
and because of the coherent nature of Josephson tun- 
neling [|T8|] . Clearly, the DJQP process allows for a far 
superior measurement of a Cooper pair box qubit than a 
SET. 
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